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1 Introduction 



In spite of a growing trend in physics to define the physical world 
as being made of information itself and thus information is de- 
fined in this way I am using this term to express some knowledge 
about things and ideas. To gain this information - before you 
need to formulate a suitable question, see Titus Lucretius Carus 
De Rerum Natura. According to Britannica Concise Encyclopae- 
dia - equations, in essence, - are questions. To get corresponding 
knowledge (information about things and ideas) we have to solve 
these equations which in many cases is not an easy task. The 
trouble is that the questions are too detailed and hence the idea 
of weak formulation of the original equations is used. In Internet 
we can find the following characterization of this idea: 

"Weak formulations are an important tool for the analysis of 
mathematical equations that permit the transfer of concepts of 
linear algebra to solve problems in other fields such as par- 
tial differential equations. In a weak formulation, an equation is 
no longer required to hold absolutely (and this is not even well de- 
fined) and has instead weak solutions only with respect to certain 
"test vectors" or "test functions"". See Wikipedia <Weak formu- 
lation>; the page last modified on 27 January 2010 by unknown 
author?. 

An example in which the idea of weak formulation is used is 
a celebrated Galerkin method. In this method the original equa- 
tions are not changed but the original spaces in which solutions are 
searched are drastically changed, for example, when the original 
space is substituted by a finite dimensional usually a low dimen- 
sional subspace. Surprising is that in this way in many cases you 
can get quite correct results even for very complicated systems 
describing, for example, the fluid flows, see, e.g.,[l], [2], [3] and 
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Internet. 

We have to remember, however, that the weak solutions of the 
weak equations, obtained in a frame of reduced-order philosophy, 
are not solutions of the original problem, see [5]. 

In the paper presented, in contrary to the canonical situation, 
we propose such use of the weak formulation of the original theory 
that the weak solutions have a clear physical interpretation and 
perhaps nice mathematical properies. So, instead of the original 
mostly partial differential equations (PDE) with the initial and 
boundary conditions (IBC) defined in a sharp way, we consider 
PDE in which IBC are trated as random or smooth quantities. As 
a consequence of that approach, the original nonlinear equa- 
tions are substituted by the linear equations for correlation 
functions or their generalizations in both cases called the n-point 
functions (n-pfs) or as in the title - the n-point information. This 
step can be treated as a weak formulation of the original nonlinear 
system of equations because considered linear equations contains 
solutions of the original theory as well as the new solutions. The 
second step in the paper is quite opposite to what is done in any 
canonical weak formulation: instead of narrowing the space 
in which n-point information are considered, we move to a larger 
space - the free Fock space (in which we do not postulate the 
permutation symmetry of functions). However, in this enlarged 
(free) Fock space - like in a smaller space of weak formulation - the 
transfer of concepts of linear algebra is possible and even general 
solutions in many cases can be constructed. 

Given the extraordinary ease of constructing unilaterally re- 
verse operations to many operators which appear in the free Fock 
space and by introducing an additional parameter (minor coupling 
constant) can be derived new equations for n-pfs. In this study 
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and others we were trying to better understand derived equations. 

In Sees 2 and 3 we define the strong and weak formulations of 
the considered equations. 

In Sec. 4 the free (super, general) Fock space metodolodgy is 
described and the basic equation for the correlation functions, 
(11) is postulated. 

Sec. 5 is devoted to the canonical perturbation theory applied 
to the Eq. 11 and a determination of the arbitrary terms which 
appear in the free Fock space, see also Sec. 8. 

It seems that very interesting is the idea of regularization of 
the true culprit of many problems of nonlinear theory: 

XiN -> Ai(I + A 2 M) i? 1 A r (1) 

whereby in the derived equations in addition to the usual sum of 
operators related to the linear ((L + G)) and nonlinear (N) parts 
of the theory, see Eq.2 and Eq.ll, the product (\2(L + G)M) 
appears, see Sees 7 and 8. In a particular case of regularization, 
M = N, - the closed equations for n-pfs - are obtained. In Sec. 9 
remarks about functional calculus are given. An example of such 
equations and the general solution for the 1-pf in the case of the 
ip s - model is discussed in Sec. 10 where some general remarks are 
also included. 

This paper is a continuation of work [9] with , as we hope, 
better use of the operator-valued functions and with useful exten- 
sion of considered operators. In comparison with work [6], where 
evolutionary type of equations were considered, in this paper - 
equations of "resovent type" are used, for comparison, see [21]. 
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2 Strong (exact) formulation 

We will assume that a theory is formulated by the following integro- 
differential equations 

L[x; <p(x)] + XN[x; ip] + G(x) = (2) 

with a linear and nonlinear dependence on the unique field (p (first 
and second terms) and a free term G. In the case of homogeneous 
environment or materials, the operators L and N do not depend 
explicitly on the vector variable x E B where a set B describes the 
domain of the unique field <p . In order to improve the description 
of equations the components of the vector x contain also discrete 
indexes which usually appear as subindices of the fields or func- 
tions. It turns out that riddance of the lower and upper indices 
is a big improvement of description. In this way only one unique 
field tp is considered. The coupling constant A, later denoted by 
Aiand called the major coupling constant, contains the memory 
of the nonlinearity of the original, strong formulation (2). This is 
usually an expansion parameter in the perturbation approach to 
the statistical and quantum fields. 
We rewrite the Eq.2 as 

(L ^)(x) + Li[x;^(x)] + AiVo[x;^(x)] + AiVi[x;^ + G(x) = (3) 

This is an equation for the unique field p{x) in which operators 
(functionals) with sub index "o" denote a local or self interaction of 
a cell or particle or field, but sub index "1" denotes an interaction 
among the constituents of the physical system. An additional 
restriction of particular terms in Eq.3 comes from a co- variant 
character of proposed equations with respect to symmetry trans- 
formations of the theory. A symmetry of equations can be used 
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to introduce another averages (or smoothing) than the ensemble 
averages with a simple, geometrical interpretation and often di- 
rect measured. It is interesting that these two kind of smoothing 
procedure lead to identical equations for n-pfs, [?],[?]• This is a 
happy coincidence because in this way we can compare theory 
with experiment without resolving the ergodic problem, [?]. The 
problem of calculation of the time averages is related to solving 
equations upon n-pfs with appropriate additional conditions. 

3 Weak formulation can be linear, multitime and non- 
commuting. Positivity conditions 

These three features of presented here weak formulation are cho- 
sen not for a provocative purpose or to illustrate a philosophical 
doctrine that science is a matter of convention but to draw the at- 
tention of practically oriented reader that they together can also 
be applied in numerical methods. Linearity is associated with 
randomization or/and smoothing of description, multitime is also 
related to a more complete randomization description in which 
the time is not distinguish and is treated as other variables. In re- 
sult, the Kraichnan-Lewis multitime correlation functions instead 
of the one time Reynolds' or Hopf 's correlation functions are used: 

y{x) ^< if(xi) ■ ■ ■ cp(x n ) >; n = 1, 2, oo (4) 

Finally, the noncommuting variables are associated with an 
additional generalization of the arena in which systems are de- 
scribed: We introduce the free Fock space in which the corre- 
lation functions or smoothed n-point functions (n-pfs) need not 
be permutation symmetric. To stress this fact we will call n-pfs 
the n-point information. In this last step we do not project the 
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original equations, in our case equations for n-point information, 
but we enlarge the space in which solutions are searched. This 
is exactly opposite to what it is done in the Galerkin methods. 
Nevertheless, like in the Galerkin methods, where corresponding 
space is diminished(l), this permits the transfer of concepts of lin- 
ear algebra to solve considered equations. Among these concepts 
we take the generators for Cuntz algebra, right and left invertible 
operators and plenty projectors using of which allows us to con- 
struct varies final formulas for generating vectors generating the 
correlation functions or, in the general case, n-pfs, which we call 
the n-point information. In the case of the correlation functions 
we have important restrictions: 



for n=0, 1,2... These restrictions we will call the positivity condi- 
tions, for correlation functions. The case 



means a trivial theory. 

4 The free Fock space and n-points information 

To deal with an infinite collection of correlation functions or n- 
pfs, < ip(xi) ■ ■ ■ (p(x n ) >, the one generating vector \V >can be 
introduced by means of which all these n-pfs can be reproduced. 
Using such a vector we describe the infinite system of branching 
equations for n-pfs in a compact form of one vector equation which 
can be transformed in varies equivalent and useful forms. From 
definition 




(5) 



< ip(x) 2n >=0-^<p = 



(6) 
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\v>= 

J2n=l I d ?'{n) < <P(Xl) ■ ■ ■ V(Xn) > ff(Xl) ■■ -^*(S„)|0 > +V |0 > 

(7) 

where operators fj*(x), hermitian conjugat to the operator fj(x), 
satisfy the Cuntz relations 

f)(x)ff(y) = i-5(x-y) (8) 

which mean that ranges of operators ff are pairwise orthogonal 
and in fact the expantion (7) is a generalization of the idea of 
expansion of a vector by means of a multiple orthogonal base. 
Here operator /- unit operator, 5- is a product of Kronecker's 
delta (discrete case) and vector |0> represents, using quantum 
field theory language, a "vacuum", 

r)(x)|0>=0 (9) 

see[?],[?] and [9]. Set of vectors (7) form the free linear Fock space 
F. Lack of commutation between quantities f)*in the generating 
vectors |V> (free Fock space) does not exclude the possibility that 
n-pfs < (p(x{) ■ ■ ■ tp{x n ) >are permutation symmetric, see below. 
Conditions (8) and (9) are enough to show that 

< 0|tK£i) • • • fj(y n )\V >=< cp(yi) ■ ■ ■ ^{y n ) > (10) 

Since operators fj do not commute, the above formula is able 
to retrieve from the generating vector7 also permutation non- 
symmetrical n-pfs like in the case of quantum fields (p. But a 
true reason for introducing non-commuting field fj is such that 
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operators introduced below and constructed by means of the op- 
erators 77,77* can be right or left invertible, inverses to which can 
be easily constructed. This leads to a variety of useful formulas 
for n-pfs. 

We postulate the following equations for the n-pfs < <p(xi) • ■ ■ (x n ) > 
which by means of the generating vector (7) can be described in 
a compact way: 



(L + XN + G)\V >= P \V > +\P N\V >= |0 > ph (11) 
with operators 



L = J 7} * (x)L[x; fj]dx + |0 >< 0| = 
/ 77 * (x)L(x, y)fj(y)dxdy + P Q 



(12) 



N = J r) * [z)N\z\ fj]dz + P N (13) 



and 



G = J f)*(x)G(x) (14) 

see[?],[?]. A small modification of the r.h.s. of Eq.ll is connected 
with a demand of right invertability of the operators L and N 
what force us to add terms Po\V >and APo^l^ >• For a concrete 
choice of that term, see (108). In fact, the |0 > p hj^ |0 > have to be 
used only for G 7^ and N ^ 0. It reminds us of a distant analogy 
with virtual particles of Quantum Field Theory and therefore is 
called the physical vacuum. In fact the r.h.s. of Eq.ll comes from 
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the fact that the original Eq.2 and the averaging process, <...>, 
does not say anything about zero component of the equation. 

Eq.ll means that we have chosen averages with respect to used 
additional conditions (ensemble averages). As we said in Sec. 3, in 
the case of homogeneous system, both types of averages lead to 
the same equations for the correlation functions. 

By introducing projectors P n projecting on the consecutive 
terms of the expansion (7), we can express the projection proper- 
ties of operators (12-14) as follows: 

P n L = LP n (15) 
(diagonal), where n=0,l,2,..., 

P n N=J2 p nNP m (16) 

(upper triangular), where n=0,l,2,..., see (13) and (108). 

P n G = GiVi (17) 

(lower triangular), where n=l,2,.... The operator values function 
N[z;fj] can be a polynomial or other function depending on the 
vector variable z and the operator variables fj(x) indexed by the 
vector variable x. The operator N is related to a nonlinear part 
of the strong formulation of theory (the original differential equa- 
tions (2). The operator G describes a source term with a function 
G{x) correponding to the external forces, for example. It is sym- 
tomatic that diagonal and upper triangular operators describe an 
interaction or selfinteraction of the constituents of the system and 
that lower triangular operators describe an interaction with the 
external world or quantum properties of the system (microworld) . 
The simplest diagonal operator is the unit operator 
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/ = |0 >< 0| + J fj * (x)fj(x)dx (18) 

Other diagonal operators are the projectors used in formulas (15- 
17) and constructed by means of the tensor product of vectors: 

P n = Jfj* {Xl) ••■fj* (#n)|0 >< 0\fj(x n ) ■ ■ ■ f}(xi)dX(n) 

(19) 

where Pq = |0 >< 0|. They form a complete set of orthogonal 
projectors: 

^ P ™ = aTld P ™> P n = PnSmn (20) 

71=0 

We can say that projections P n \V >, for n=l,2,..., provide n- 
points information about the local nature of the system but the 
projection Po\V >provides rather global, agregated information. 
In the case of system representing the Universe, the r.hs. of Eqs 
like (11), (59), (76), (80) can be interpreted as a vacuum, see [9]. 
Thus, in this interpretation the (classical) vacuum contains the 
global information about the Universe. Like in QFT a non-trivial 
structure of the vacuum arises only through the nonlinear theory 
and this is a positive element that might enable its stady. Identi- 
cal equations as (11) take place in QFT, for vacuum expectation 
values 

< <p{x 1 )...(p(x n ) > (21) 

where (/3are now operators with appropriate equal time commu- 
tators. In this case however, n-pfs (10) are not permutationally 
symmetric. 
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5 Right invertible linear part of the theory and approx- 
imated solutions 

If the kernel L of a diagonal operator L is a right invertible: 

J L(x, y)L R l {y, z)dz = S(x - z) (22) 

then the operator L is a right invertible in the Fock space F, see 
7-9. A right inverse to L, denoted by L i? 1 can be constructed as 

Lr = fj * {z)Lft(z, w)f){w)dzdw + Pq (23) 

where, as before, the symbol J means the sumation or integration 
with respect to components of vectors z,w. In fact, the operator 
L R satisfies the weaker equation: 

LLx 1 = I (24) 

Having constructed a right inverse to a given operator L, we 
can construct the projector on the null space of L: 

U L = I- L R l L (25) 

Multiplying Eq.ll by a right inverse operator {L + G), we can 
describe this equation in an equivalent way as follows 

[/ + X(L + cy^N] \v >= u L+G \v > + 

{L + G) R l [P Q \V > +\P N\V >) 

(26) 

where the projector on the null space of the operator (L + G) is 
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fl L+G = i-(L + G)R(L + G) (27) 

From point of view of Eq.ll, considered in the free (full, super) 
Fock space F, the projection n^+clF >can be any vector from 
space IIl + gF. With different right inverse operators!/™ 1 , different 
projectors Ul + gou the null space of the operator (L + G) are 
constructed. In result, a general solution to the generating vector 
|V> obtained by means Eq.26 has different parametrizations. 

From definition of smoothing operations given in Sec. 3, we see 
that obtained n-pfs are permutation symmetric. Denoting by S a 
projector on vectors generating permutation symmetric n-pfs , we 
should have, for the physical solutions: 

\V >= S\V > (28) 
Hence and from Eq.26 , we get 

[I^XS(L + G)^N]\V>= 
SI1 L+G \V > +S(L + G)^ 1 (Pq\V > +\P N\V >) 

(29) 

Eq.29 was derived by using first Eq.28 and next acting on Eq.26 
with projector S. Similar procedure is used in the Galerkin method, 
however, there is an important difference with Galerkin approach, 
namely - Eq.28 - is an exact equation and hence the solutions to 
Eq.29 can also be exact solutions of the original problem. 

Taking into account the projection properties of operators S, L^ 1 
and G, Eq.ll can be also described as: 

{/ + A(J + SL^Gy^L^NjlV >= 
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Sfl L \V > +S{L + G) R l (P \V > +\P N\V >) 

(30) 

This or Eq.29 are the output equatios for calculating succes- 
sive approximations to the generating vector |V> expanded in the 
positive powers of the coupling A standing at the nonlinear part 
of original theory (2): 

a 

\v >=J2 XJ \ V > {]) ( 31 ) 

3=0 

The arbitary element of Eq.26, TIl\V ^lF, can also be ex- 
panded in this way 

oc 

^l\V>=Y, XJ (^l\V>P (32) 

3=0 

The terms of the above series can be restricted by the permutation 
symmetry of n-pfs. In the case of symmetrical solutions (28) to 
Eq.29 we can assume that the symmetry part of the arbitrary 
element Iix|V >is given by a linear theory: 

Sfl L+G \V >= S\V > (0) = \V > (0) (33) 

It is a common situation accompaning to almost every approxi- 
mated and exact theory see also Sec. 8. It means that only due 
to terms with A ^ higher approximations appear. So the zeroth 
order approximation 

\V > (0) = SU L \V > (0) (34) 
The first order approximation 
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\V > (1) = -{(I + SL^G^SL^NjlV > (0) (35) 
The second order approximation 

\V > (2) = -{(/ + SL^G^SL^N^V >W (36) 

and so on. All these approximations can be obtained from a single 
vector formula: 



\V >= 

A A A A -| A /V A -| A A A A -t A A A 

{/ + A(/ + SL R l G)- l SL R l N}-\I + SL R l G)- l SU L \V > 

(37) 

and it is not inconceivable that these compact formula provides a 
new look at old divergent problems of unrenormalizable theories. 

Let us focuse on the zeroth order approximation, (34). It is a 
symmetrical solution to the Eq.ll, for A = 0. We will consider this 
equation with additional simplification that exterial forces acting 
on the system and represented by the operator G = 0. So we get 
the following vector equation 

L\V > (0) = (38) 

describing a linear theory. In the "components" form it looks as 
follows: 

J dx x L{x, xi) < ip(xi) ■ ■ ■ (f{x n ) > (0) = (39) 

for n=l,2,... Let us assume that we know a function A which 
satisfies equation 
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dxiL(x,xi)A(xi,y) = 



(40) 



for all x,y. With functions A, we can construct, for even n, 
symmetrical solutions to Eqs 39 as follows: 



and zero, for odd n. Because of symmetry, the permutation oper- 
ation is taken in the sum. In a graph representation, with vertexes 
denoted by xi, x n: a particular terms of sum 41are represented 
by graphs with one edge vertexes. It is easy to see that, for such 
solutions, conditions (5) are satisfied. In addition, we also expect 
that these conditions are also fulfilled for the higher order ap- 
proximations for the correlation functions when the perturbation 
parameter, the coupling constant A, is a small quantity in Eq.30. 
As an example of 41, 



< <p(x)<p(y) > W= -(A(x, y) + A(y, x)) = A(x, y) (43) 



where the permutation symmetry of function A is assumed. Hence 
we see that in zeroth order approximation the correlation functions 
(41) are sums of products of the 2-point function43. 

To find a connection of the zeroth order 2-pf A with the general 
solution to Eq.2, in the case of A = 0, G = 0, we represent this 
solution in the form 



<i P (x 1 )---i P (x n )>^=V(x 1 ,...,x n )^ = 

h. Y.perm A (^l> ^) A(x 3 , X±) • • ■ A(x n _l, X n ) 



(41) 




(42) 



1(3 



^«] = /^(W) (44) 

in which the kernel T of the general solution (2) (A = 0, G = 0) 
represents possible elementary solutions to Eqs (2) and40 - labeled 
by the vector variable u which are summed with u- dependent fac- 
tor a. We do not define here a set to which the vector parameteru 
u belongs. From 43) and 44 we get 



< <p(x)(p(y) >(°)= A(x,y) = 
J dudw5aT(x, u)T(y, w)a(u)a(w)P[a] 

Denoting the functional integral occurring in (45) as 



(45) 



/ 5aa(u)a(w)P[a] = P(u,w) (46) 
we get for the function A, the following expression: 

A(x,y) =< (f(x)ip(y) >^= J dudwT(x,u)T(y,w)P(u,w) 

(47) 

which relates the correlation function A between cells with the 
kernel F of the general solution to (2), at A = 0,G=0, and the 
second order moments P(u,w) of the probability density P[a\. 

The kernel T of the general solution to Eq.2, in case (A = 
0, G = 0), can be constructed by means of a right inverse L^ 1 , 
see 22: 

T(x, u) = Hl(x, u) = 5(x — u) — J dzL^(x, z)L(z, u) (48) 
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In such case this is a projector (idempotent): 

and the subspace described by the projector T can be identified 
with all solutions to Eq.2 . In many cases, with every subspce, 
one can define a linear differential equation, see [13]. In fact a 
differential equation can be represented geometrically by the set 
of all solutions, and a symmetry of the differential equation is 
defined as a map which transforms this set into itself, [14]. 
For a symmetrical T 

T{x,y) = T{y,x) (50) 

and the Gaussian 



P(u, w) = 5(u — w) (51) 

we get from (45) 

A(5, y) =< <p(x)<p(y) >(°>= T(x, y) (52) 

In other words, for the linear systems (2) described by the symmet- 
rical, idempotent kerneirof the general solution and the Gaussian 
type moments (51), the zeroth order 2-pf is identical with this 
kernel (sic!). It is astonishing case which shows that, at least at 
the zeroth order level, information lost by averaging or smoothing 
procedures (< <p(x) >^= for all x) can be recovered with the 
help of 2-point correlation function (52). It is not excluded that 
the following hypothesis is true: - in the case of nonlinear equa- 
tions a similar phenomen takes place with the help of the higher 
order correlation functions and the Gauss smoothing. 
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For (^nonlinear, local theory, the above properties are par- 
ticularly important because they allow to substitute the bilinear 
products A • Aoccurring in such a theory by the one A, see Sec. 5. 

It is interesting to notice that condition (50) called sometimes 
the reverse normalized condition is in fact one additional con- 
dition imposed on a right inverse operator to get a solution to 
(2) with smallest Euclidean norm (generalized inverse or pseudo- 
inverse called also Moore- Penrose inverse). This illustrates the 
kind of processes described by Eq.52. It turns out that every 
bounded operator L : K — > H (Hilbert spaces) with closed range 
has a generalized inverse. To see properties of such operators, 
see[15]. 

6 Right invertible "nonlinear" part of the theory 

In this case we will assume and it actually takes place in the free 
Fock space, [6]- [9] that a right inverse operator exists to the 
operator N: 



and that the operator can be explicitly constructed in the 
most general form, [6]- [9]. Introducing a projector on the null 
space of the operator N: 



NN P = I 



(53) 



P = I — N R l N = 1- Q N - < 0|7V = 



(54) 



we can write equivalently Eq. 11 as follows: 



[/ + X^Nr^L + G)]\V >= P\V > + 




(55) 
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With symmetry (28) we get 



[I + X^SN^iL + GW >= SP\V > + 

X^SN^ 1 (P \V > +XP N\V >) (56) 

Assuming that with limes A — > we get finite results for the 
generating vector |V>, we have 

\V(X = oo) >= SP\V(oo) > +F iv|^(oo) > (57) 

In a spirit of perturbation theory we assume that the arbitrary 
term of solutions to Eq.56 

SP\V >= SP\V(oo) >= \V(X = oo) > -P N\V{oo) >=? 

(58) 

and a possible answer to that is a real issue. 
7 Modified or regularized theory 

We mean by this an introducing to the term representing the 
nonlinear part of the original theory - one extra parameter: N — > 
N(X2). This parameter is called the minor coupling constant and 
is denoted by A2 where the original parameter A, now denoted by 
Ai, will be called the major coupling constant. If at the end of 
calculations we go with parameter A2 to zero and with A2 — > 0, 
N(X2) — > N we will speak about a regularization of the theory in 
other case about its modification. In certain cases such theories 
with two coupling constants Ai, A2 7^ satisfies the Dyson crite- 
rion of convergency,[16], and hence the name - regularization. In 
other words, divergences of the usual perturbation theory can be 
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interpreted as residues of the deformed, nonlocal and more ba- 
sic theory which is currently still unknown. For other interesting 
ideas concerning renormalizability of a theory see [17]. 

So, instead of Eq.ll we consider the resolvent regularized or 
modified equation 



(L + Ai(J + \ 2 M)rN + 6) \V >= 

P \V > +A 1 P (/ + \ 2 M) R l N\V >= |0 >reg (59) 

where the operator (I + X 2 M) R 1 is a right inverse operator to the 
operator (7 + AM). In the mathematical literature, this operator 
is called a generalized resolvent of the operator M , [11]. See [21], 
for other use of expression - the resolvent equations. 

The r.h.s. of Eq.59 comes from the fact that the original Eq.2 
and the averaging process, <...>, does not say anything about zero 
component of the equation. Therefore, acting on this equation 
with projector Po we should get identity. 

Multiplying the last equation by the operator (I + AM), we 
remove, at least from the l.h.s. of Eq.59, ambiguities contained 
into the right inverse operator and get 



[I + \ 2 M)(L + G) + \iN 



\V > 



I + A 2 M) ( P \V > +X 1 P (I + \ 2 M) R l N\V > 



(60) 



This equation can also be treated as a regularization of the original 
Eq.ll in a sense that, for \ 2 — > 0, (60) tends to Eq.ll. 

Eq.60 is not equivalent to (59} because it was derived by multi- 
plying by the right invertible operator (I — X 2 M). This operation 
is equivalent to the operation of multiplication by the projector 
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Q(i+x 2 M) = (I + MM) R \i + A 2 M) (61) 

We have here an interesting situation in which the projected Eq.60 
tends in the limes, A 2 — > 0, to the original Eq. 11. Is it possible 
from the point of view of logic? The answer to this question is 
positive, if we take into account that equation (4.5) considered in 
the Fock space are overdetermined. Why? Because Eqsll or its 
regularized version (59) or (60) are exactly the same in the free 
Fock space as well as in the physical space (28). 

Now we transform the Eq.60 in an equivalent way taking into 
account the right invertibility of the operators M and L. We get 

{ J + [A 2 M(L + G)]- R l [L + G + AiiV] } \V >= P M \V > + 
[A 2 M(L + G)] R l (P \V > +AiP (/ + X 2 M) R 1 N\V >) 

(62) 

where 

P M = I- [M(L + G)]- R l M(L + G) (63) 

is a projector on the null space of the operator M(L + G). This 
projector has the following structure: 

Pm = I-(L + G) r 1 Qm{L + G) (64) 
with projector Qm = M R l M and a right inverse 

{L + G) R l = L R \i + GL R l ) R l (65) 

The projected vector of |V>, Pm|^ >, is an arbitrary vector of 
the null space of the operator M(L+G). It describes an additional 



22 



freedom which a deformed theory brings to the original theory 
given by Eq. 11. From (27) we have 



Pm&l+g = Ul+g (66) 

It is also possible to have a different perspective on the Eq.59. 
This time we will not weaken this equation but we represent the 
operator R(X 2 ) = (I + X 2 M) R 1 in a more explicit way assuming 
that M is a right invertible operator. Then we get 

R(X 2 ) = (J + X 2 M) R l = 
[I + (A 2 M) i? 1 ]- 1 [(A 2 M) i? 1 + f M R(X 2 )} (67) 

with an arbitrary element f mR(X 2 ) where the projector 

f M = I - M R l M = I-Q M (68) 
It is interesting that a choice 

f M R(\ 2 ) = (69) 

leads to a generalized resolvent R(X 2 ) satisfying the resolvent 
equation as 

■R = -MR 2 (70) 



d\ 2 

which is also satisfied for the operators M with the usual resol- 
vents. This equation, in turn, ensures that, for at least right 
invertible operators M, the choice R=[I+ (X 2 M) R 1 ]- 1 (X 2 M) R 1 
satisfies the same condition as for the non singular M, namely: 
R — > 0, for A 2 — > 00. In practical considerations, the term 
TmR{X 2 ) 7^ can be used when there are problems with limes 
A 2 ^0. 
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Let us also notice that in the zero- \2limit the vector 



|tt(A 2 ) >= [/ + (A 2 M)fl 1 ]- 1 (A 2 M)fl 1 |$ >^ |$ > (71) 
at the assumption 

A 2 M|*(A 2 ) >^ 0, for A 2 -> (72) 
because from (71) results that 

[A 2 / + M R 1 ]\^(X 2 ) >= M^|$ > (73) 

and finally 

[A 2 M + /]|*(A 2 ) >= |$ > (74) 
In other words the operator 

[i+{\ 2 M) R l ]-\\ 2 M) R l ^i (75) 

behaves as the unit operator, with A 2 — > 0, if the vector M|^(A 2 ) >has 
oscillations suppressed by A 2 — > 0. In fact, the paper presented is 
about the indeterminate expression • 00 or rather -0 • do- because 
we are dealing here with operators. 

Eq.59 with representation (67) is following: 

(l + g) \v > + 

Ai ([/ + (X 2 M) R 1 ]- 1 [(X 2 M) R 1 + f M £(A 2 )]) N\V > 

= P \V > +X 1 P Q (I + X 2 M)- R l N\V > (76) 

With a right inverse operator, (L + G)^ 1 , we can transform this 
equation in an equivalent way as 
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|l/>+A 1 (L + G) i? 1 [/ + (A 2 M) i? 1 ]- 1 - 
([(A 2 M) i? 1 + f, /J R(A 2 )])iV|y > 

= tl L+G \V > + 
(L + G)^ 1 (P \V > +A 1 P (/ + \ 2 M) R l N\V >) (77) 

where a projector on the null space of the operator (L + G) is 
given by 

U L+G = i-(L + G) R 1 (L + G) (78) 

With condition (28), Eq.77 leads to an equation in which the 
permutation symmetry of solutions is explicitly taken into ac- 
count: 

\V > +X l S{L + G) R 1 ■ 
([/ + {\ 2 M) R l ]- l [{\ 2 M) R l + f M R{\ 2 )]) N\V > 

= Sfi L+G \V > + 
S(L + GQ" 1 (P |^ > +XiP (I + \ 2 M) R l N\V >) (79) 

For a choice of (95), this equation in which the arbitrary el- 
ement - the projection 11^+^1^ >- is the same as in the theory 
with major coupling constant Ai = 0, see Sec. 5, can be closed, see 
the next section. It is noteworthy that perhaps equations: (62), 
(79) are a new type of equations which offer a new perspective 
on the calculation of n-pfs, see Sec. 10. It is encouraging that in 
both formulations unspecified elements can be chosen in such a 
way that similar results can be obtained. 
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There is possibility of a different regularization of Eq.ll: 

[L + AiiV(/ + \ 2 M)- R l + 6) \V >= 
P \V > +X 1 P N(I + X 2 M)~ R l \V > (80) 
Hence and from (67) we get 

{l + +g}|f> + 

AiiV ([/ + (X 2 M) R 1 ]- 1 [(X 2 M) R 1 + f M R(X 2 )]) \V > 

= P \V > +X 1 P Q N(I + X 2 M)r 1 \V > (81) 
With choice (95) 

\V > +(L + • 
(&[J + (AsiV)" 1 ]- 1 + X 1 N ([/ + (X 2 N)- 1 }-'f N R(X 2 )})) ■ 

■\v>= ii L+G \v > + 

(L + G)^ 1 (p |V > +X 1 P Q N(I + A 2 iV) i? 1 |y >) (82) 

With assumption (28) the above equation leads to its symmetrized 
version 

\V >+5(L + G)^- 
(&[J + (A^ 1 ]- 1 + A^ ([/ + (A 2 iV)- 1 ]- 1 f i v J R(A 2 )])) • 

•|y >= sn L+G |v> + 

S(L + G)^ 1 (p |V > +A!P iV(i + A 2 ivVl^ >) 

(83) 
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This equation as other previous equations, for the choice (95), 
leads to closed equations for n-pfs. A positive fact is that, as in 
perturbation theory, the undetermined vector STI^+dV >can be 
identified with the original linear theory (N = 0) and that for 
derivation of the last equation we did not project, besides their 
symmetrization, the previous equations. This last fact makes that 
weakening of the original Eqs (2) still can be interpreted as the 
appropriate averaging. 

The arbitrary term, TnR(\2) 7 can be chosen in such a way to 
remove the "secular" terms which arise when the A2 — > 0. 

Finally, we can connect two kind of regularizations and consider 
the following regularization: 

N -+ \ (N(I + A 2 iVV + (J + X2N) r 1 n) (84) 

to mimic in some degree a commutation of N with operator (/ + 
A2./V) if the last operator would exist. 

For such regularization and choice (69) chosen only for simplic- 
ity, we get the following equation: 

/ + tA L + g )r (t 7 + (^r^IQn + i})}- 

■\V >= SI1 L+G \V > +S(L + G) R l ■ 
(P \V > +±P (N(I + A 2 iV) i? 1 + (/ + \2N) R 1 N^) \V >) 

(85) 

We see that at the diagonal terms the rate of coupling constants 
Ai^appears and particularly, the 1-pf V only has such a depen- 
dence. This means that, for monomial theories at least and used 
assumptions, the 1-pfs V, in the \im\i — > 0, stop to depend on 
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the major coupling constant Ai. If we agree that exact closure of 
equations for n-pfs V is equivalent to a summation of the infinite 
terms of the canonical perturbation theory, we can compare this 
result with disappearing of the effective coupling constant of the 
perturbation theory. This rather negative result can be treated as 
a hint that A2 7^ 0, and that non polynomial theories should be 
taken into consideration, or, that Aidepends on A2in such a way 
that y const when A2 — > (some kind of renormalization of 
the original theory). We can imagine the following situation in 
which the developed type of theories can be used: For a certain 
class of initial and/or boundary conditions ("laminar" conditions), 
phenomena are well described by the polynomial theory (2) with 
A2 = 0. For others ("turbulent" conditions, we need to use aver- 
ages or smearing and then non polynomial terms (A2 7^ 0) have 
to be taken into account because, as we see in averaged version 
of Eqs (2), Eqs (59) and (60), A2enters these equations in the in- 
verse powers. In other words, the non polynomial terms become 
apparent only when considered situations (boundary and initial 
conditions) are becoming increasingly complex and averaged so- 
lutions and correlation functions are needed. 

Let us show that at least for a choice (69) the right invertible 
operator M dose not commute with its general resolvent R = 
[/+(AM)^ 1 ]- 1 (AM)^ 1 : We get 



M, R 



A- 1 / A- 1 



jM = 



= A _1 f 



M 



I+(\M)r 

see (64) and (68), for definitions of appropriate projectors. 
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8 Undetermined terms of the general solution. 



If we confine ourselves to the solutions satisfying the condition 
(28), then the projected by the projector S, Eq.62, can be written 
as follows: 

{/ + S[X 2 M(L + G)] R l [L + G + AiM \V >= SP M \V > + 
S[X 2 M(L + G)] R l (P Q \V > +A 1 P (/ + \2M) R l N\V >) 

(87) 

Introducing the operator 

A = I + S[X 2 M(L + G)] R l [L + G] (88) 
we rewrite Eq.87 as follows: 

{A + Ai5[A 2 M(L + G)]^N}\V >= SP M \V > + 
S[X 2 M(L + G)]^ 1 (P |V > +AiP (/ + X 2 M) R l N\V >) 

(89) 

This equation may serve as a starting point for a perturbation the- 
ory of the modified Eq59with the major coupling constant Aias an 
expansion parameter and the operator S[X 2 M (L + G)] R N as the 
perturbation operator. With assumption about non singularity of 
the operator A, we can rewrite Eq.?? as 

{/ + Aii- 1 S , [A 2 M(L + G)]^n\ \V >= A- l SP M \V > + 

A- l S[X 2 M{L + G)] R l (P \V > +A 1 P (i + X 2 M) R l N\V >) 

(90) 
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We will assume that, for Ai = 0, the theory is continuous and 
hence, the unknown term 

A- l SP M \V >(°) +A- l S[X 2 M{L + G)] R l P Q \V >(°) = 

\V]\ 1 = Q>=\V >^ (91) 

where the generating vector \ V > (^generates correlation functions 
of the originally linear theory characterized by the zero coupling 
constant Ai. In a general case , for Ai 3= 0, the undetermined 
element of Eq.90 is chosen in a such way that the IBC are satisfied. 
Additionally, if the expansion (31) is used and we get 

A~ l SP M \V\ Ai, A 2 > +A- l S[\ 2 M{L + G)]^P \V > = 

\V > (0) + E°li K(A- l SP M \V; A 2 >)0 (92) 

then we try to choose the undetermined projections SPm\V; A 2 >^ 
in such a way that no divergent terms do not appear in the se- 
ries (31), when A 2 — > 0. We have assumed here that the original 
theory (A 2 = 0) is well defined. 

The simplest case is then again when 

A-^PmI^ > +A- 1 S[X 2 M(L + G)]^P |V > 
= \V > (0) o 

&P M |F > +S[X 2 M(L + G)]^P |^ >= i(A 2 )|V > (0) 

(93) 

This assumption is in agreement with (91) and also corresponds 
to the case when the IBC are satisfied by the zero-order approxi- 
mation and no divergences appear in the expansion with respect 
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to the major coupling constant Ai. The assumption that arbi- 
trary term of Eq.90, A~ 1 SPm\V >, is identified with the original 
linear theory, see Eq.2, indicates that the further consequences of 
this theory are related to the nonlinear effects (perturbation of 
the original, linear theory, (Eq.2)) occurring in the system. Quite 
natural assumption which is a started point of almost every per- 
turbation theory. Again we must recall that assumption (93) does 
not exclude that the part of arbitrary element Pm\V >G PmF 

(i-S)P M \V>=f(\ h \ 2 ) (94) 

responsible for the elimination of asymmetric parts in Eq.62 can 
be a complicated vector function allowing to satisfy the condition 
(28). 

If at the choice (93) some bad, "secular" terms appear in the 
power series expansion with respect to the major coupling con- 
stant Ai, see [18] then to remove such terms we can use the for- 
mula (92). It appears an open question: Is it possible using this 
type of equations to gain a new perspective on the renormalization 
in of quantum and statistical field theories? 

Next we assume 

M = N (95) 

because such a choice of the modified operator often leads to closed 
or partly closed equations for n-point correlation functions. It is 
natural that a candidate for modified theory should be looked 
among closed theories of the type (95). It is worth noting that a 
more general choice of the modifying operator 

M = AN (96) 
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with a diagonal operator A can also be used to close the considered 
equations. We assume that operator Ais a non singular operator 
not to increase freedom of the theory contained in the projection 
TmR(^2), see (67). Nevertheless, the choice of (95) exhibits a 
kind of self-similarities of fractal theories. 

9 Functional calculus and considered equations 

The functional calculus is define sometimes as a branch of math- 
ematics about inserting operators into functions to get in result 
meaningful new operators, see, e.g., [19], [20]. Such operator- 
valued functions often appear in solutions to many linear equa- 
tions of physics and engineering problems and in our opinion are 
an important generalization of the concept of function. 

From the functional calculus point of view we can say that our 
paper needs a definition of the function (/ + \ 2 M)~ in the case 
of upper triangular operators M like (95). The natural and often 
used definition of this function as the power series 

(I + \ 2 M)- 1 := I - X 2 M + (A 2 M) 2 - ... (97) 

for unbounded operators, in every its term can be incorrect and, 
moreover, in the case of upper triangular operators M very incon- 
venient because it introduces to every correlation function of basic 
Eq59 an infinite number of modification or regularization terms. 
Therefore, we used a different definition 

(I + X2M)- 1 := (J + X 2 M) R 1 (98) 

where subscript "R" means a right inverse operator to the oper- 
ator / + A2M. So, instead of a senseless expression (97) we use 
well defined and explicitly constructed expression (98), see (67), 
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(95) and (107)-(109). The operator (98) is called the generalized 
resolvent of the operator M . Interesting thing is that regulariza- 
tion or modification of equations (11), given by the generalized 
resolvent of the operator M = N, leads to closed equations for 
the correlation functions. Moreover, many operator-valued func- 
tions in the free Fock space can be defined by means of the type 
of Dunford-Cauchy integral with generalized resolvent (98): 

/a 2 (M) = / /a 2 (A)(I + \M) R l d\ (99) 
r 

It is not unlikely that this paper is the first step to understand- 
ing these operator valued functions in the context of equations for 
the correlation functions. We can say that we analyzed such a 
problem: to which conclusions leads a regularization or modifi- 
cation of Eq.ll which consists in replacing an upper triangular 
operator N by the diagonal + lower triangular operator, e.g., 



N -+ f X2 (N)N = [ I f X2 (X)[I + (XNftl^XNtfdX ) N 

(100) 



or 



N^\(fx 2 (N)N + Nfx 2 (N)) (101) 

see (67), (69) and (95). In the last formula like in (84) we took 
into account that different right inverse operators do not commute 
with each other: Denoting two different right inverses by N R and 
N R if we assume their commutativity we come to the absurd: 
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N([N R \ N R r m >= N R V \^ > -N^l® >= (102) 

for any vector |$ >. Also the operator valued functions con- 
structed by different right inverses do not commute with each 
other and the original operator N. For example, we expect that, 
at least formally, iVcommute with operator valued functions like 
(/ + AiV) _1 or exp(-N) and so on. Nevertheless, we try to mimic 
commutativity of their formal predecessors (formal inverse oper- 
ators, see (97)). We can support this idea in the following way: 
in the case in which the problem in consideration is ill-posed and 
we use averages, then the original Eq.2 looses its meaning and a 
more free connection between averages and Eq.2 is recommended. 
In this way we substitute traditional theories with very compli- 
cated perturbation series and related to the open equations for 
correlation functions by a new, more flexible theories with closed 
equations like (??). 

In the case of regularization of the original theory with nonlin- 
earity described by the operator N we have a demand 

limx 2 ^of\ 2 (N) = I (103) 
or its weaker form, e.g., 

Um x ^ofx 2 (N) = / - P (104) 

It is not clear whether the condition (103) or (104), although 
weaker than the condition (71) and (72), can be fulfilled for the 
operator valued functions occurring in formulas (100) and (101). 
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10 An example of the theory with a local nonlinear 
term, divergences and final remarks 

Let us consider the equation: 

J L(x,y)if(y)dy + Xif 3 (x) = (105) 

which is called the (p^- or Hurst model. In this model the function 
G = in Eq.2. Usually, the kernel L is a sum of the Dirac's deltas 
and its derivatives or their discrete version. For this model, the 
operator 

L = y f) * (x)L(x,y)f)(y)dxdy + Po (106) 

the operator 

N = j \*{z)f){~z) 2 dz + PqN (107) 

where, to guarantee right invertibility of the operator N, we have 
chosen 

P Q N = PoNPi =Poj dzr\{z)dzPi ^ (108) 
In this case a right inverse to the operator N can be chosen as 

n r 1 = Jifj * (y)) 2 v(y)dy + Nr'Po (109) 

with 

N^Po = PiNrPo = Pif)*(y)Po (HO) 
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We will use the regularized theory leading to closed equations, 
Sees 7-8. For a sake of simplicity, we use Eq.79 in the case of 
(95) and the operator G = 0. We also choose (69), for the arbi- 
trary element of the generalized resolvent (67). In result, we get 
equation 

{ J + X 1 S(L)- R 1 {[I + (X 2 N)- R 1 ]-\X 2 N)- R 1 ) N} \V >= 

^n^i^ > +S(L) R l (P \V > +X 1 P (I + X 2 N) R 1 N\V >) 

(111) 

and further 

[i + ±S(L) R l {[I + (X 2 N) R 1 ]- 1 Q N ) } \v >= Sfl L \V > 
+S'(L) i? 1 (P \V > +X 1 P (i + X 2 N) R 1 N\V >) 

(112) 

where projector 

Q N = N R 1 N (113) 

If projection STI^T/ >can be calculated from Eq.112 in which 
Ai = we would get 

✓s A A A -t A / r\ \ A A /r\\ 

sn L \v > ={/-SL fl 1 Po}|v> <0 >={/-p }|v><°> 

(114) 

From Eq.112 we get then 
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s\ s\ \ A A -| f A -t A \ A 

P 2 sii L \v >= A|v > (0) (ii5) 

and 

s\ s\ /\ 1 A A -I A A A A A / f\\ 

Pi{I + ^SL^Q N }\V > = P^U^V >= P^V >(°) 

(116) 

where we took into account that the operator N R is lower trian- 
gular and 

PoQn = (117) 

These are closed equations for the lowest - 1 and 2-pfs. We 
can't get an equation for the 0-pf, Vq, see (7). With Eq.116 we 
can express the 1-point solution as: 

P 1 \V >= A{/ + y'SL^Qn}- 1 ^ > (0) (118) 

To describe explicitly these equations we have to take into ac- 
count that 

Q N = N-'N = 
J r)*(x) 2 fj{x) 2 dx + P 1 J ff(z)pQfi(y)dyPi (119) 

and 

Lr-Qn = irf(x)L^(x 1 y)ff(y)f)(y) 2 dxdy + 
P 1 L] i 1 P 1 ff)*(z)P r}(y)dyP 1 

(120) 
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with arbitrary value for the variable z. From (116) we get the 
following equation for the 1-pf: 



V(x) + ^ J L- R \x,y)dyV(z) = V^°\x) (121) 
with arbitrary fixed value for the variable z . Denoting the integral 



L R \x,y)dy = a(x) (122) 



we rewrite Eq.121 as 



V(x) + ^a(x)V(z) = V^(x) (123) 



Its solution is 



V(x) ee V(x; z) = V^\x) - \(x) V ^1 {Z 1 (124) 

The result is rather frustrating because, for A2 — > 0, dependence 
of the 1-pf V on Ai disappears completely. Similar results we get 
for the 2-pf with trivial condition (6), for A2 — > 0. This could 
mean that perhaps we should keep A2 7^ , or, we should try to 
use other regularization also proposed in Sec. 7. 
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